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A simple model of Dark Matter (DM) which couples with the two nucleons with different coupling 
strengths without violating isospin conservation is explored as an example of the importance of 
keeping the experimental data analysis as model independent as possible. The apparent contrast 
between the recent CRESST event excess and current XENON100 spin-independent (SI) exclusions 
is shown to be removed by simply not assuming that the couplings of the DM with the two nucleons 
are equal. 



The identification of Dark Matter is a long standing 
problem of modern astroparticle physics and cosmology. 
A particle candidate to be a constituent of DM needs 
to be stable on the time scale of the Universe age, have 
a suitably small annihilation rate at its typical densities, 
and the strength of its interaction with photons and bary- 
onic matter must be weak enough to explain the absence 
of clear traces of its presence other than the gravitational 
effects on astronomical or cosmological scales, like the 
marked separation of the gravitational lensing from the 
x-ray and light emission in cluster collisions ( More- 
over, since a high number density would make the annihi- 
lation rates high, the DM particle(s) should be massive, 
reasonably heavier than the proton. Since no particle 
in the Standard Model of Particle Physics combines all 
these characteristics, the particle theory of DM has to 
be an extension of the Standard Model. Experimentally, 
the search is carried on either indirectly by looking for 
secondaries of the DM annihilations (like 7 or positron 
excesses and vs) , or directly, by seeking rare nuclear scat- 
tering events in underground detectors. Heuristically, the 
lowest order effective scattering interaction has a Spin- 
Independent (SI) component given, in the nonrelativistic 
limit, by the product of the number densities, and a Spin- 
Dependent (SD) component arising, if the DM particle 
has spin, from the dot product of the spins of the DM 
particles and the target nuclei. 

But, while for the SD component it is generally recog- 
nized that the coupling strengths with the two nucleons 
may be different, for the SI component it has become 
customary to assume, in the analysis of direct DM ex- 
perimental results, that the couplings f v , n with the two 
nucleons are equal. This assumption is a priori arbitrary, 
restricting the analysis to some models, like the MSSM 
lightest neutralino with squarks much heavier than the 
higgs (if the latter condition is not satisfied, the MSSM 
neutralino interacts with f p 7^ /„ [2]). In fact, a general 
lowest order effective lagrangian for the SI scattering of 
a DM fermion \ on " nucleons is 0, 0] : 

£ <* XXifpPP + fnnn) (1) 
whose scalar operators xXi PPi an d nn act, in the non- 



relativistic limit, as particle number density operators. 
The resulting spin-independent zero-momentum transfer 
cross section for the elastic scattering of a DM particle 
off a nucleus is then determined by the particle numbers: 



a^fi 2 (f p Z + f n (A~N)) 2 



(2) 



The experimental reason for reporting limits extracted 
from a null result assuming f p = f n is that this scenario 
is the best constrained, or, equivalently, the most easily 
detectable, by any direct DM experiment, since the scat- 
tering cross section of Eq. ([2} is maximal for any target 
nucleus. This allows to compare the discovery potential 
of experiments which only reported limits. But, when 
comparing limits to a positive signal (or two positive sig- 
nals), which may arise from a more generic scenario, f p 
and /„ should be left independent of each other. This is 
known to reduce or even remove the apparent incompati- 
bility of results like CoGeNT and XENON @,@], provided 
that suitable values off p , n are employed. Moreover, it is 
possible to constrain [4|], and even determine, f p and f n 
by combining different experiments. 

A simple example of a DM interaction producing dif- 
ferent scattering cross sections on the two nucleons is the 
following minimal extension of the standard model, that 
exhibits f p =/= f n while manifestly conserving isospin. The 
new (dark) particle content consists of an isospin spinor 
doublet = (ip, x), along with a scalar isospin triplet 
<f> and a scalar singlet p mediating the dark interaction. 
The lagrangian is: 

C = g t $ ■ + NaN) + g s p{^ + NN) (3) 

where a are the Pauli matrices for the isospin vectors, 
and g t ^ s are the triplet and singlet coupling constants. 
N = (p, n) is the nucleon doublet. Such a scenario could 
arise from a weakened strong interaction in the dark sec- 
tor, resulting in highly massive dark mesons coupled to 
both the nucleons and the dark fermions with the weak- 
ened g s ,t constants. If M s ,t are the respective masses 
of the scalar singlet and triplet, the low energy effective 
lagrangian for elastic scattering is: 
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c = [rjj-j ■ NSN + [wj ^ NN ( 4 ) 

If we now assume that the dark doublet mass splitting 
— M x is large relative to the DM kinetic energies, 
with the lightest dark spinor being x (e.g. — M x ^> 
1 - 10 MeV, for 1 < M x < 1000 GeV), then the inelastic 
scattering is suppressed, reducing Eq (g]) to: 



= - (jfc) XXNa 3 N + (§A XX{nn+pp) = 
\Wt) XX&P - nn) + (jj-^ XX(nn + pp) = ( 5 ) 
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which clearly gives Eq © with f p = g 2 /M 2 - g 2 /M 2 
and /„ = g 2 /M 2 + g 2 /Mf , and assumes the lightest dark 
spinor constituting the DM is \- Of course, if the lightest 
dark spinor has positive isospin (ip), f p = g 2 /M 2 +g 2 /M 2 
,-md/,. ,/: M 2 gf Mf. 

The terms gt<f> ■ ^<J^ and g s p^^> in Eq. Q provide 
annihilation channels of the dark fermions ^ to the p and 
(j> scalars, as well as all standard model particles lighter 
than x or i 3 to which p and <fi may couple. Moreover, 
whenever M tjS ^ M^ iX , the effective lagrangian for an- 
nihilation to any isospin doublet A has the same form of 

Eq. (gl), i.e. (-§^j ^ct* ■ Act A + (jfc) W/VA, which 
has enough free parameters to easily provide the correct 
DM relic density in a way analogous to what shown in 
Ref. 0- E.g., if the lightest dark fermion is x, the annihi- 
lation cross section of x to an isospin singlet pi (xX ~^ W 4 ) 
in the center of mass reference frame is: 



g 2 M x (M x - Mft) _ g 2 M 2 



Mf 



Mt 



(6) 



where the last expression is valid for M n <§; M x . As for 
the impact of the model-blind SI analysis with f p ^ /„ on 
current direct dark matter results, it has been shown in 
the literature 0,0] that by simply choosing f p /f n «-(Xe 
neutrons) /(Xe protons) ~ — 78/54 « —1.4 it is possible to 
suppress the zero momentum transfer cross section of the 
Xe isotopes enough to make the exclusions from XENON 
compatible with the CoGeNT positive signal, and, at the 
same time, reconcile the 2 annually modulated positive 
signals currently available. This, of course, applies to 
whatever model can yield essentially arbitrary / Pi „, in- 
dependently on whether they violate isospin conserva- 
tion/symmetry or not. It is therefore possible that the 
DM has an isospin dependent interaction (f p ^ f n ) with 
the nucleons without even violating isospin conservation, 
and these scenarios can better fit the known observa- 
tions and exclusions than the more theoretically biased 



fp = fn models. Indeed, at present, the most appropriate 
approach to DM data analysis and interpretation is still 
to simply drop the f p — f n restriction, waiting for the 
more recent positive signals to become solid (CoGeNT 
admittedly still has to perform many systematic studies, 
and CRESST 0] doesn't yet have a modulation) to ex- 
tract from the data clues on which theoretical model best 
fits the observations. 

As a new example of how f p ^ /„ can reconcile pos- 
itive signals and exclusions (provided that the detector 
targets are different), and of how the data from differ- 
ent experiments can tell us which the possible values of 
/ Pi „ are, let us consider the latest DM signal claimed by 
CRESST-II, focusing on the M2 candidate of Ref. [|, 
which, in the standard f p — f n analysis, has M x = 11.6 
GeV and a = 3.7 x 10 -2 fb, in clear contrast with both 
the w 3.3 x 10~ 4 fb limit of the XENON100 and the 
w 5.2 x 10~ 3 fb of the CDMS-II exclusions Now, 
if an experiment reports an f p = f n cross section <jn 
for the nucleon, it means its isotope (Z,A) had a zero 
momentum transfer DM cross section 



u A = UNA 2 [ t± ) 2 
Pp 



(7) 



where pa,p are the DM-nucleus and DM-proton reduced 
masses, and the small difference between the 2 nucleon 
masses has been neglected. But from Eq. ([2]) it is clear 
that 



_ Al 

ffN [Z+f n /f p (A-ZW 
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(8) 



f p /f n Z+{A-Z)* 



where, of course, ct„ is the neutron, and <t/v the nucleon 
cross section. Picking the most abundant isotopes of each 
elements in the respective detectors, and using the above 
values of ctat, it is easy to obtain from Eq. approxi- 
mate proton and neutron cross sections for XENON100 
and CRESST-II in the scenario f n /f p = -0.7 of Ref.s 



'p 



16 fb for 132 Xe of XENON100, a p w 4.2 fb 
for 74 Ge of CDMS-II, a p « 1.6 fb for both 16 and 40 Ca 
of CRESST-II. Clearly, in this Xe sensitivity suppress- 
ing scenario and this rough approximation, the signal of 
CRESST-II is not excluded by XENON100 nor CDMS- 
II. In fact, also the neutron interaction is not excluded: 
ct„ « 8 fb for Xe 132 of XENON100, ct„ w 0.8 fb for both 
O 16 and Ca 40 of CRESST-II, and the a n « 2 fb for 74 Ge 
of CDMS-II is still compatible with the 11.6 GeV can- 
didate of CRESST-II, though the maximal suppression 
for 74 Ge occurs at f n /f P = —0.76, equal and opposite to 
the Z/(A — Z) ratio of 74 Gc. This is why, in this rough 
estimate and specific scenario, the CDMS-II limit seems 
closer to the CRESST measurement than XENONlOO's. 
The reason for the identity of ct PjJI of 16 and 40 Ca 
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CRESST-II intersection with XENON100 and CDMS-II allowed ellipses CRESST-II intersection with XENON100 and CDMS-II allowed ellipses 




FIG. 1: CRESST-II (blue thick line), XENON100 (violet horizontally hatched area) and CDMS-II (red hatched area) for 
M x = 11.6 GeV (left), and detail of the upper left intersection region of the 3 experiments (right). As the overlap region of 
the red ellipsis allowed by CDMS-II and the violet ellipsis allowed by XENON100 intersects the CRESST-II elliptical contour, 
there are values of f p and f n for which the 2 results are compatible. Not surprisingly (see text), these values have a ratio fn/ fp 
close to -0.7 (dash-dotted line) which maximally suppresses the SI sensitivity of the Xe isotopes. 



in CRESST-II is that, for all isotopes with A = 2Z, 
A 2 /[f P /f n Z + {A- Z)] 2 = A 2 /[Z + f n /f p (A - Z)Y = 
4/(l + /„// P ) 2 . 

To take a more rigorous approach, I'll follow Ref. 
with some modifications. The ellipses displayed in Fig. 
Q]are given by Eq.s (6,7) of Ref. [4|], with the following 
differences: 

- instead of the coupling strengths, the adopted vari- 
ables are the signed square roots of cr p , n , so that 
the sign ambiguity in the first of Eq.s (7) (Ref. [1]) 
can be removed, while the right hand term remains 
1. 

- For CRESST-II, the equation is an equality, while 
for XENON100 and CDMS-II it remains an in- 
equality (limit). This means that XENON100 and 
CDMS-II allow the interior of the respective el- 
lipses, while CRESST-II affirms the proton and 



neutron cross sections for M y 



11.6 GeV are, 



where dR/dE is the rate per unit target mass, a a the 
scattering cross section of isotope A, pa the DM particle- 
nucleus reduced mass, FF\ the nuclear form factor, Ja 
the abundance of isotope A in the target's composition 
and I a = / f(v)/vdv the quasi moment -1 of the DM ve- 
locity distribution relative to the detector, which depends 
on A via the lower integration extremum (minimum in- 
cident velocity able to produce the recoil energy E). The 
rightmost side of Eq. ([9]) employs the simplifying as- 
sumptions of neglecting the weak dependence of I a on A 
{I A ~ I) and that for light DM the momentum transfer 
is small (so that FF\ « 1). From Eq. ([2]) we can deduce: 



a a [ZJop + (A - Z) y /on]' 1 



= A 



2 (7jV 

a 2 
I'p 



(10) 



within errors, somewhere along its ellipsis. 



again containing the approximation /i p rs fj, n , the right- 
most side being the value of the center one for the 
signed JcTp~ = ^/a^, i.e. for f p = f n . Substituting Eq. 



(pTO]) in the rightmost side of Eq. ((9J and dividing by 
p/2M x IJ2 A IaA 2 ^, we obtain: 



The conversion from <7jv to a p ^ n starts from the stan- 
dard rate equation: 



2J 



/ap-Z + ^iA- Z)f 



Ea/aA 2 



= 1 



(11) 



dR _ p ^ OA 2 , r 



(9) 



2M V ~ 4^ P? A 



which produces the ellipses of Fig. [T] The consistency 
of the above approximations has been verified by inter- 
secting the ellipses with the Ja^ = ^JaZ. line, finding 
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CRESST-II excluded by XENON100 but still allowed by CDMS-II 
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FIG. 2: CRESST-II (blue thick line), XENON100 (violet solid 
area) and CDMS-II (red hatched area) for M x = 25.3 GeV. 
As XENONlOO's ellipsis is entirely contained in CRESST-II's, 
the Ml candidate is incompatible with XENONlOO's exclu- 
sions at this DM mass for any SI couplings. CDMS-II is 
instead still compatible with CRESST-II for f p /f n 1.2. 



that the coordinates of the intersections are indeed the 
square roots of the <t/v reported by the experiments for 
M x = 11.6 GeV. 

XENON100 and CDMS-II allow, for M x = 11.6 GeV, 
the inside of the respective ellipses, while CDMS-II claims 
the most likely values of ^/Op and yja~/l are along its el- 
liptical contour. Since the latter intersects the overlap 
area of the XENON100 and CDMS-II allowed regions, 
there are some values of ^/dp and y^ovT from CRESST- 
II still allowed by both XENON100 and CDMS-II. The 
left Fig. [T] shows the two symmetric areas compati- 
ble with the CRESST-II M2 candidate, while the right- 
hand plot shows one of them in detail. Given that the 
fn/ fp = —0.7 crosses the middle of the two symmetric ar- 
eas still allowed, the /„/ f p ratio is constrained to values 
near ~ —0.7. 

For completeness, let us consider also CRESST-II's 
candidate Ml (M x = 25.3 GeV, « 1.6 x 10~ 3 fb, vs 
the fa 1.2 x 10~ 5 fb limit of XENON100 and w 2 x 10~ 4 
fb of CDMS-II in the f p = f n scenario), which has a 
69% contribution from tungsten. From the latter cir- 
cumstance, it is reasonable to expect that XENON100 
constrains CRESST-II more strongly than at 11.6 GeV, 
since the (A — Z)/Z ratios of tungsten's isotopes are 
closer to those of xenon than they are for oxygen and 
calcium, so that scenarios suppressing XENONlOO's sen- 



sitivity to the Ml candidate would largely suppress also 
CRESST-II's sensitivity. Fig. [2] confirms this expecta- 
tion: XENONlOO's exclusions are incompatible, at max- 
imum likelyhood, with the Ml candidate, for every choice 
of fp/fn- CDMS-II, instead, is still compatible, for 

fp/fn 1-2, with Ml. Exploring the full 2a Ml region 

of Ref. [8| is beyond the scope of this paper, and would 
probably require the detailed knowledge of CRESST-II, 
XENON100 and CDMS-II data that only the respective 
collaborations have. 

In conclusion, it is still easy to devise simple models of 
DM where the DM particle couples with the two nucleons 
with different coupling strengths, even without violating 
isospin conservation, and can yield the correct cosmolog- 
ical relic abundance of DM. This circumstance makes it 
necessary, for a correct comparison of experimental re- 
sults some of which are positive signals, to analyze (or 
reanalyze) the data without assuming f p = f n . In fact, 
while for comparing the achieved sensitivities of experi- 
ments without positive signal it makes sense to pick the 
scenario in which all experiments are most sensitive, pos- 
itive signals may well arise from a scenario where this is 
not the case. The latest example of the above is the 
apparent contrast between the recent CRESST-II event 
excess and current exclusions, like those of XENON100 
and CDMS-II, which are currently compared in a model 
biased fashion. Dropping the f p = f n bias shows that 
only the Ml SI candidate is excluded for all f p . n , while 
for M2 f p /f n is constrained to values close to -1.4. 
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